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In the current cosmological model, only the three lightest elements were cre- 
ated in the first few minutes after the Big Bang; all other elements were pro- 
duced later in stars. To date, however, heavy elements have been observed 
in all astrophysical environments. We report the detection of two gas clouds 
with no discernible elements heavier than hydrogen. These systems exhibit the 
lowest heavy-element abundance in the early universe and thus are potential 
fuel for the most metal poor halo stars. The detection of deuterium in one 
system at the level predicted by primordial nucleosynthesis provides a direct 
confirmation of the standard cosmological model. The composition of these 
clouds further implies that the transport of heavy elements from galaxies to 
their surroundings is highly inhomogeneous. 

In modem cosmological theory, the light elements and their isotopes are produced during 
the first few minutes after the Big Bang when the universe cools during expansion from temper- 
atures T ~ 10^ K to below ~ 4 x 10*^ K. In this brief epoch, termed Big Bang Nucleosynthesis 
(BBN), D, '^He, ^He, and ^Li are synthesized with an abundance ratio relative to hydrogen that 
is sensitive to the cosmic density of ordinary matter (i.e. the baryon density ^b,o)- BBN theory 
also predicts negligible production of the heavy elements with abundance ratios X/H < 10^^° 
and one must await the physical conditions that are typical of the stellar interiors (i). 

The analysis of gas observed in absorption along the lines-of-sight to high-redshift quasars, 
distant galaxies that host supermassive black holes, is a powerful probe of the BBN yields. Par- 
ticular attention has been given to deuterium, partly due to observational convenience but also 
because the D/H abundance ratio is very sensitive to ilb,o- For quasar sight lines, the measured 
log(D/H) = -4.55 ± 0.03 (2, 3) translates into r];,,o/i^(BBN) = 0.0213 ± 0.0010, which is 
fully consistent with the value inferred from the Cosmic Microwave Background (CMB) power 
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spectrum Qb^oh^iCMB) = 0.02249lg:gggp (WMAP7 (4)). This excellent agreement between 
two essentially independent experiments stands as a marked triumph of the Big Bang theory. 

On the other hand, all of the systems with measured D have heavy-element abundances 
that exceed, by many orders of magnitude, the BBN prediction. In fact, despite measurement 
of thousands of galaxies from the early universe (e.g. (5, 6)) and careful study of the diffuse 
gas that permeates the universe (e.g. (7, 8)), one has yet to detect anything near primordial 
enrichment. For structures denser than the mean cosmic density (Fig. 1), the high-redshift 
universe has exhibited a floor in the metallicity Z, the mass fraction of elements heavier than 
hydrogen, at ~ 1/1000 of the solar abundance (Zq). Similarly, among several old and iron-poor 
stars, only one has a metallicity Z ~ 10"^^© (9), with all the remaining having enhanced C or 
O abundances. The existence of a minimum level of enrichment at about 1/1000 solar has been 
associated with the metal production in Population III (PopIII) stars, primordial stars that form 
in metal free environments via H2 cooling. In fact, models and numerical simulations (e.g. 
(10, 11)) show that ejecta from this first stellar population can enrich the interstellar medium 
(ISM) of the host halos up to ~ 10~^Zq and pollute the surrounding intergalactic medium 
(IGM) as soon as 1 billion years after the Big Bang {z ~ 6). 

In this work, we report on the hydrogen and metal properties of two gas 'clouds' at z ~ 3, 
when the universe was only two billion years old. We observed quasars SDSS Jl 13418.96-1-574204.6 
(zem = 3.522) and Q0956 + 122 (^em = 3.297) on UT January 3 and 5, 2006 and UT April 7, 
2006 with the HIRES spectrometer on the Keck I telescope on Mauna Kea (see SOM text 1). 
Previous low-resolution spectra of these quasars from the Sloan Digital Sky Survey had shown 
significant absorption at wavelengths A < 4000A characteristic of substantial optical depth tll 
at the H I Lyman limit (at wavelength A < 912A), typical of the Lyman limit systems (LLSs). 

Our Keck/HIRES spectrum of Jl 134-1-5742 reveals a sharp break in the flux at ~ 4000 A, 
indicating the presence of a LLS at z ~ 3.4 with tll > 2 (hereafter named LLS1134; Fig. 
2), confirming the lower resolution data. A search for absorption lines at a redshift consistent 
with the Lyman break further reveals the presence of the H I Lyman series through to Lyman- 
22 (H I 913.5), corresponding to A ~ 4030A. In the high resolution spectrum, two distinct 
absorbers can be identified within LLS 11 34, the main system (hereafter LLS 11 34a) dX z — 
3.410883 and a weaker component (named LLS1134b) dX z — 3.41167, separated by 5v ~ 
54 km s~^. In contrast, two flux decrements in the normalized HIRES spectrum of Q0956 -|- 122 
(Fig. 3) are visible at ~ 3860 A and ~ 3750 A, revealing the presence of two LLSs with tll ~ 1 
at ^ ~ 3.22 (LLS0956A) and 2 ~ 3.10 (LLS0956B). 

A closer inspection of these spectra also reveals no detectable metal-line absorption at the 
velocity of the strongest H I component for neither LLS 1 134a nor LLS0956B. Such a complete 
absence of heavy element absorption has not been previously reported for data with comparable 
sensitivity. For LLS 1 134a, besides the H I Lyman series, the only other detected transitions are 
D I Lya and Ly/3, offset by the appropriate 5v = —82 km s^^ from the hydrogen absorption. 

To characterize the physical properties of these LLSs (summarized in Table 1), we modeled 
the hydrogen and deuterium absorptions using a minimization algorithm (Table SI; Fig. S2 
and S3, see SOM text 2) and we derived upper limits on the column densities of various ioniza- 
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tion states of heavy elements (Table S2; Fig. S4, see SOM text 3). To translate these limits into 
constraints on the gas metallicity, it is necessary to make assumptions on the ionization state of 
the gas (see SOM text 4). Every LLS with A^hi ^ lO^^cm"^ analyzed to date has exhibited 
absorption characteristic of a predominantly ionized gas (e.g. {12, 13)). Standard interpretation 
is that the medium has been photoionized by an external radiation field, presumably a combi- 
nation of the extragalactic UV background (EUVB) generated by the cosmological population 
of quasars and galaxies together with emission from local sources (e.g. a nearby galaxy). This 
conclusion is based on comparison of the observed ionic column densities of heavy elements 
with simple photoionization models. Every previous LLS has shown substantial absorption 
from doubly and triply ionized species, e.g. 81+''" and C+^, which trace predominantly ionized 
gas. Parametrizing the ionization state in terms of the ionization parameter U = $/(c rin), with 
$ the flux of ionizing photons and nn the gas volume density, all previously analyzed LLSs 
have exhibited U > lO^'^ corresponding to < 10~^ cm~^ at 2; ~ 3 (Fig. S6). 

Adopting this U value as a limit to the ionization state and current estimates for the spec- 
tral shape of the EUVB {14), we infer metaUicities of Z < IQ-^-^Zq and Z < IQ-^-^Zq for 
LLS 11 34a and LLS0956B (see Fig. S5). These upper limits are 100 to 1000 times lower than 
typical measurements of LLSs and over an order of magnitude lower than any previous observed 
metallicity at 2; > 2 (Fig. 1). These limits are only comparable to the abundances detected in 
the most metal poor star (9) and are suggestive of a primordial composition. Super-solar metal- 
licity is commonly found in the surroundings of quasars as well as in a few LLSs. MetaUicities 
between solar and ~ 1/10 solar are typical of galaxies, and sub-solar metal enrichment down to 
IO^^Zq is characteristic of the ISM and IGM at early epochs. Remarkably, the most iron poor 
stars in the Galactic halo, thought to be the repository of the first generation of metals, have 
total heavy element abundance comparable to these limits, but generally above ~ 1/1000 solar. 
Our analysis uncovers regions of the universe at 2; < 6 with essentially primordial enrichment, 
whose traces can be found within the oldest stellar populations in the present universe (9). 

The detection of deuterium in the metal free LLS 1 134a provides a direct confirmation of the 
BBN. From the analysis of the D and H absorption lines, we derive log(D/H) = —4.69 ± 0.13 
(see SOM text 2). The observed value is consistent at la with the theoretical predicted value 
log(D/H) = —4.592 {15), assuming from WMAP7. This measurement is in agreement 
with previous determinations in quasar absorption line systems at z > 2 (Table S3; Fig. 4). 
Because the lack of metals in LLS1134a confirms its pristine composition, this agreement 
strengthens the hypothesis that at low metaUicities (IO^^Zq or less) the observed deuterium 
abundances are representative of the primordial value {16) and astration cannot be responsi- 
ble for the lingering scatter in the observed D/H. Deuterium abundances from quasar ab- 
sorption line systems are therefore solid anchor points for models of galactic chemical evo- 
lution (e.g. {17)). Combining D/H in LLS1134a with values from the known D-bearing sys- 
tems, we obtain a logarithmic weighted mean log(D/H) = —4.556 ± 0.034 that translates into 
fife 0^^ (BBN) = 0.0213 ± 0.0012, after accounting for both random and systematic errors {3). 
Consistent with previous studies, we do not include in the weighted mean the error on the as- 
sumed level of the quasar continuum light. 
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The absence of metals in LLSl 134a and LLS0956B is outstanding also in the framework of 
theories for the metal enrichment of cosmic structures (Fig. 5). Numerical simulations suggest 
that LLSs typically arise in galaxies (e.g. (18)), in dense gas above ~ 100 times the mean 
cosmic density Pmean- At the same time, in models of the IGM enrichment, metals are ejected 
to hundreds of kpc from star forming regions, resulting in substantial pollution of the nearby 
gas. As a consequence, the metallicity predicted for stars or for the ISM at z < A (19) ranges 
between 0.1 — IZq, three orders of magnitude higher than the limits inferred for these two LLSs. 
Similarly, the metallicity predicted for the hot halo of galaxies and for the surrounding IGM 
(19, 20) exceeds by a factor of 10 or more the limits for LLSl 134a and LLS0956B. Contrary to 
any prediction and any previous observation, these two LLSs reside at significant overdensity 
(p/ Pmean < 850) but in an unpolluted portion of the universe. 

Metallicity below 10~^Zq appears even exceptional when compared to the level of pre- 
enrichment from PopIII stars that is predicted by models (10, 21, 22) already 500 million years 
after the Big Bang (z ~ 10). Our limits place additional constraints on the widespread dispersal 
of metals from primordial stellar populations and the first generations of galaxies. Remark- 
ably, the gas we detected could in principle fuel PopIII star-formation at z ~ 3, because its 
metallicity lies at or even below the minimum enrichment required for metal cooling to induce 
fragmentation in the collapsing material in the absence of dust (23). Given a lower limit on the 
LLS physical size £ = A^hi/ (xmnn), with xni the neutral fraction and the total hydrogen 
volume density, we can infer the total hydrogen mass in these clouds Mh = mpiVn^^. We find 
Mh > 6.2 X 10^ M0 for LLSl 134a and Mr > 4.2 x 10^ M© for LLS0956B, comparable to the 
mass of the mini-halos where the first generation of stars formed (e.g. (22)). Therefore, if this 
gas were able to collapse further and shield from the ambient UV radiation, it would potentially 
give rise to PopIII stars two billion years after the transition between the first and second gener- 
ation of stars (PopII) is thought to have occurred. Thus, pair-production supemovae associated 
with the death of these massive and metal-free stars may be found even at modest redshifts (24). 

The pristine composition of LLSl 134a and LLS0956B can be reconciled with model pre- 
dictions and previous observations if mixing of metals within the IGM is an inefficient and 
inhomogeneous process. A varying degree of metal enrichment is seen in multiple components 
of LLSs (12, 13, 25), implying that mixing does not operate effectively on small scales. The de- 
tection of ionized metals in LLSl 134b, the weaker component at +54 km from LLSl 134a, 
reinforces this point. Further, studies of metal systems in the low density and diffuse IGM sug- 
gest that at least some of the ionized heavy elements (e.g. C IV) are in small and short-lived 
clumps (26). A low volume filling factor for metals is also consistent with theories of metal 
ejection from supemovae, in which most of the heavy elements are initially confined in small 
bubbles (27) and only subsequently diffuse in the surrounding IGM. Plausibly, LLS 11 34a and 
LLS0956B originate in a filament of the cosmic web where primordial regions coexist with en- 
riched pockets of gas. If the metal enrichment is highly inhomogeneous, these two LLSs could 
just be the tip of the iceberg of a much larger population of unpolluted absorbers that trace a 
large fraction of the dense IGM. 

Beside the implications for the BBN and the metal distribution, the detection of metal free 
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LLSs is tantalizing in the context of galaxy formation and evolution. Modem theory and sim- 
ulations predict that most of the gas that sustains star formation is accreted in galaxies through 
dense and narrow streams, known as cold flows (28,29). These gaseous filaments are highly ion- 
ized by both the EUVB and the radiation escaping from the central star forming regions. Cold 
flows should therefore appear as LLSs in the spectra of bright quasars (30-32). Further, in con- 
trast to metal enriched gas that is outflowing from galaxies, this infalling material is expected to 
be metal poor (31). Although direct observational evidence of cold flows is still lacking, these 
streams are thought to be ubiquitous at high-redshift and primordial LLSs such as LLS1134a 
and LLS0956B are ideal candidates for this elusive mode of accretion. 
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LLS1134a 



LLS()956B 



Redshift 

log A^HI 
logD/H 



3.410883 ± 0.000004 
17.95 ±0.05 
-4.69 ±0.13'^ 
15.4 ±0.3 
< (1.43 ±0.05) X 10^ 



3.096221 ± 0.000009 
17.18 ±0.04 



6Hi(kms ^) 
Temperature (K) 
Metallicity (Zq) 



< 10-4-2 

< -2.10 

< -1.86 
> -3 



20.2 ±0.8 
< (2.48 ±0.19) X 10^ 



< 10-3-8 

< -2.40 

< -1.98 
> -3 



logxm 

lognn 

logC/^ 



" Including LLS1134a, the current best estimate for the primordial deuterium abundance becomes log(D/H) = 
—4.556 ± 0.034. ^ The Usted values are physically motivated, but not directly measured. Note that the metalUcity, 
ojhi, and rin depend on the assumed value. 

Table 1: Summary of the physical properties for LLS1134a and LLS0956B. For each system, we 
present: redshift, the hydrogen column density (Nm), the deuterium abundance (D/H), the Doppler 
parameter (6hi), the temperature, the metalUcity, the hydrogen neutral fraction (xm), the total hydrogen 
volume density (nn), and the ionization parameter (U). 
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Figure 1: Peak to peak variation of the observed metallicity in z > 2 cosmic structures at different 
densities. The blue, green and red rectangles show the spread in observed metallicities for diffuse gas in 
the universe; respectively, these are the IGM, LLSs, and galactic ISM. Orange rectangle: H II regions 
in galaxies. Light blue rectangle: quasar broad line regions. The black point with error bars marks the 
mean metallicity and the standard deviation for z > 1.5 LLSs. Galactic halo stars (stars connected with 
a dashed line) are represented at arbitrary density for visualization purposes. The top axis translates the 
number density in the overdensity above the mean baryon cosmic density at z = 3.5. The upper limits 
on the metallicity for LLS 1134a and LLS0956B are shown with green arrows, assuming log{7 = —3. 
Higher ionization parameters would shift these limits to lower densities and lower metallicity. See the 
SOM text 5 for additional details on the observations presented in this figure. 
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Figure 2: Top panel: Keck/LRIS spectrum of the QSO J1134+5742. A LLS at z ~ 3.4 is clearly visible 
from the break at ~ 4000 A. Middle panels: H I Lyman series transitions in the LLS 11 34a. Superim- 
posed to the data are the best-fit model (red line) and the 2a errors (yellow shaded regions). Individual 
components included in the model are marked with thin gray lines and the position of the hydrogen 
and deuterium are indicated by vertical dotted lines. The second hydrogen component (LLS 1134b) at 
+54 km s^^ is marked with a dash-dotted line. Bottom panels: selected strong metal-line transitions in 
LLS 1134a (dotted lines) and LLS 1134b (dash-dotted lines). Unrelated absorption from the IGM con- 
taminates the Si III transition (see Fig. S4). 
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Figure 3: Top panel: Keck/HIRES spectrum of the QSO Q0956 + 122. The flux decrements at ~ 3860 
A and ~ 3750 A reveal two partial LLSs at z ~ 3.22 (LLS0956A) and z ~ 3.10 (LLS0956B). The 
corresponding Lyman limits are marked with vertical lines. Middle panels: hydrogen Lyman series 
transitions for LLS0956B shown relative to z = 3.096221. Superimposed to the data, the best-fit model 
(red lines) and 2a uncertainties (yellow shaded regions). Bottom panels: selected strong metal-line 
transitions for LLS0956B. 
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Figure 4: Deuterium abundances as a function of metallicity ([X/H] = log(X/H) — log(X/H)0) for the 
z > 2 absorption line systems. Green symbols are for LLSs, while red symbols are for higher column 
density absorbers (A'^hi > 10^*^ cm^^). Metallicities obtained with silicon are indicated by squares, 
while those obtained with oxygen are indicated with circles. The right-hand axis translates the deuterium 
abundance to the cosmic baryon density JI^ q^^- The inferred fife (BBN) is shown with a solid blue 
line, together with the la errors (light blue shaded area). The CMB value and \a errors from WMAP7 
are instead shown with a solid black hne and a gray dashed area. 
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Figure 5: Overview of analytic models and simulations for the metal enrichment of the universe. Light 
blue and green shaded regions: IGM metallicity from PopIII stars with mixing between 1 — 0.05 and 
different star formation histories (10,21). Red dashed region: critical metallicity that marks the transition 
between PopIII and PopII stars (23, 33). Orange line: analytic model for the gas metal content in the 
universe from PopII stars and galactic winds (20). Gray shaded region: gas metallicity within halos 
from hydrodynamical simulations that include yields from both PopIII and PopII stars (22). Blue shaded 
region: analytic model for the ISM metallicity at different halo masses (10^^ — lO^^M©) and different 
wind models (34). Black lines: metallicity from hydrodynamical simulations with momentum driven 
winds (19) in condensed gas (solid line), hot halo (dotted line), warm-hot intergalactic medium (dashed 
line), diffuse gas (dash-dotted line), ISM (dash-triple-dotted line), and stars (long-dashed line). Upper 
limits on the metalhcities of LLSl 134a and LLS0956B are marked with green arrows. 
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Supporting Online Material 



1 Observations and Data Reduction 

Previous spectra of quasars SDSS Jl 13418.96+574204.6 and Q0956 + 122 had shown substan- 
tial absorption at wavelengths A < 4000A characteristic of the continuum opacity of H I gas. 
LLS 1134a and LLS0956B were identified as part of an ongoing program to explore the inci- 
dence and physical characteristics of LLSs (72, 35) from a sample of ~ 50 absorbers, many of 
which have apparently low metallicity (< 1/100 solar). 

We first observed quasar Jl 134-1-5742 on UT January 5, 2006 with the HIRES spectrometer 
(36) on the Keck I telescope on Mauna Kea. The data were obtained in two exposures totaling 
6300 seconds, with the instrument configured with the red cross-disperser, the C5 decker, and 
the kv380 blocking filter to suppress second order flux. In this configuration, the instrument 
provides a resolution of ~ 8 km s^^ FWHM, and the echelle and cross-disperser angles were 
set to provide wavelength coverage 4017 A < A < 8534 A. These spectra have gaps in coverage 
at A ~ 5450A and 7050A due to spacings between the CCDs in the detector mosaic. There are 
also a series of gaps in wavelength coverage beyond A > 6300A where the free spectral range 
of the spectrometer exceeds the detector width. 

We observed Q0956 + 122 with HIRES once on UT January 3, 2006 for 7200 seconds, and 
again on UT April 7, 2006 for 1800 seconds. For the January data, HIRES was configured with 
the UV cross-disperser and the C5 decker, again providing ~ 8 km s^^ FWHM resolution. The 
wavelengths covered in this configuration were 3621 A < A < 5533 A, with gaps in coverage 
only related to the CCD mosaic. The April 2006 data had HIRES configured with the red cross- 
disperse, the CI Decker, and the og530 blocking filter. The wavelengths 5465 A < A < 9998 
A were covered with wavelength gaps due to the spacings of the CCD mosaic and the limited 
detector size. 

For both quasars, the 2D images were reduced with the HIRESredux pipelineQ which ex- 
tracts and coadds the data. An interactive continuum fitting procedure within the pipeline as- 
signs a continuum level to each coadded spectral order, and the continuum normalized spectra 
are then combined into a single ID spectrum. For Jl 134-1-5742, the resultant ID spectrum has 
a signal-to-noise ratio S/N ~ 25 per 2.6 km s~^ pixel at 6300A, the central wavelength. For 
Q0956 + 122, the ID spectra have S/N ~ 30 per 2.6 km s'^ pixel at 4500A for the UV cross- 
disperser data, and ~ 10 per 1.4 km s^^ pixel at 6000 A for the red cross-disperser data. 

We also observed Jl 134-1-5742 on UT 2011 July 1 with the LRIS spectrograph (37) on 
the Keck I telescope on Mauna Kea. We obtained a single exposure of 400 seconds, with the 
instrument configured with the 600/4000 grism, the D560 dichroic and the 600/7500 grating, 
tilted to ensure continuous spectral coverage between the blue and red arms. The 2D images 
were reduced with the LowRedux pipelinea which calibrates and extracts the data. 

^ http ://w w w.ucolick. org/^xavier/HIRedux/index. html 
^http://www.ucolick.org/~^xavier/LowRedux/index.html 
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2 Hydrogen and Deuterium Analysis 



We now describe the H I and D I analysis for each LLS. 

2.1 The Lyman Limit System at z = 3.410883 toward J1134+5742 

Both the strong Lyman Umit absorption visible in the Keck/LRIS spectrum (Fig. 2), and the 
saturation of the high-order Lyman series lines in the Keck/HIRES data (Fig. S2) imply a large 
H I column density, A'hi. Fig. SI shows a zoom-in of the Keck/LRIS spectrum for Jl 134-1-5742, 
focusing on the Lyman limit absorption. We estimate the flux just redward of the Lyman limit 
to be 13 ± 2 in our relative units. This estimate appropriately includes the attenuation of the 
quasar light by the Lya forest. In the following, we assume a flat spectrum for the continuum 
and note that a more realistic spectrum would yield a slightly higher estimate for the A^hi value. 
Just blueward of the Lyman limit (A = 3950 — 4000A), we set a conservative upper limit to 
the relative flux of 0.05. Combining these two measurements, we estimate a lower limit to the 
optical depth at the Lyman limit of tll > ln(9./0.05) > 5 for a total A^hi > 10^''''^ cm"^. 
The model for this A^hi limit, shown as the green line in Fig. SI, is marginally acceptable at 
A ~ 4000A, but predicts a recovery of the flux at A < 3800A that appears inconsistent with 
the observations. Metal-lines and strong Lya lines in the HIRES spectrum, however, reveal 
the presence of other absorbers along the sightline to Jl 134+5742 at redshifts z — 3.0283, 
z = 3.0753, z = 3.1047, z = 3.2396, z = 3.2616, and z = 3.3223. These systems could 
contribute to the opacity at A < 3942 A. If they do not, then LLS 1 134a must have a higher A^hi 
value. 

Fig. SI shows the data binned with a weighted mean in windows of 50A between 3400 — 
4000A. Overplotted on these binned evaluations (whose errors only reflect statistical uncer- 
tainty) is a model curve for A^hi = io^^ °5 cm~^ (blue line). This A^hi value reproduces the 
observed absorption to 3400A. Larger A^hi values are not ruled out by the analysis, espe- 
cially when one allows for systematic error (e.g. sky subtraction). In summary, we report a 
lower limit to the total A^hi value of 10^^^ cm^ and note the data require yet higher values if 
LLS 11 34a dominates the opacity at A < 3700A. 

To further refine the A^hi measurement and its redshift distribution, we have performed 
detailed analysis of the Lyman series lines. This analysis, however, is challenged by the fact that 
gas with A^Hi ~ 10^® cm^^ has Lyman series lines that lie on the saturated portion of the curve- 
of-growth (e.g. (12)). As such, we proceed cautiously as this may impose large uncertainties in 
the model. We may first set a strict upper limit to the total A^hi of LLS 1 134a from the absence 
of substantial damping wings in the Lya profile (Fig. 2). Because the flux would be over- 
absorbed for all models with A^hi exceeding lO^^^cm"^, this is the largest A^hi value allowed 
for LLS 1134a. Second, the complete absorption at A ~ 4030 A by the Lyman series requires 
A^Hi > lO^'''^ cm~^, consistent with the Lyman limit analysis performed above. 

Lastly, we can model the H I absorption using the VPFIT v9.5 package that allows one to fit 
multiple Voigt profiles to the spectrum. Our general approach is to include as few components as 
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possible (each defined by a redshift z, column density A^hi, and Doppler parameter 6hi) to match 
the Lyman series absorption. The data are best described by two principal H I components 
(LLS1134a and LLS1134b), as evident from a visual inspection of the higher order Lyman 
series lines (Fig. 2 and S2). Since LLSl 134b is unsaturated beyond Ly-14, we first measure the 
line-parameters of this subsystem using the Lyman series lines between Ly-10 and Ly-16. We 
find z = 3.41167 ± 0.00001, 6hi = 18.0 ± 0.9 km s"^ and log A^hi= 16.71 ± 0.02 indicating 
that this component has a nearly negligible contribution to the observed Lyman limit. Adding 
an additional uncertainty of 0.02 dex related to the 1 — 3% error in the continuum placement, our 
best fit value for the column density of LLS 1 134b becomes log iVHi= 16.71 ± 0.03. Then, we 
model LLSl 134a and the deuterium absorption, using both high order Lyman series lines (Ly-7 
through Ly-16) as well as Ly-a and Ly-^. For these transitions, we must include absorption 
from gas unrelated to LLS 11 34 which we ascribe to coincident Lya absorption from lower 
redshift absorbers. A summary of the model parameters is given in Table SI. 

The strongest H I component has z = 3.410883 ± 0.000004, log ArHi= 17.94 ± 0.05 and 
an H I Doppler parameter 6hi = 15.4 ± 0.3 km s~^. The latter is consistent with a predomi- 
nantly ionized gas at a temperature T < 1.4 x 10^ K. The line-parameters for LLS 11 34a are 
well constrained by the full Lyman series and the resultant A^hi value is consistent with our 
Lyman limit analysis of the Keck/LRIS spectrum. Despite the line-saturation, there is substan- 
tial constraint on the 6hi value, and thereby the A^hi value, from the line-profile shapes of the 
Lyman series lines. This conclusion hinges, however, on the assumption that the H I absorption 
at 2; ~ 3.41088 is dominated by a single component. As we will demonstrate in the next para- 
graph, this assumption is well supported by the analysis of neighboring D I absorption. From 
the above discussion, we conclude A^hi = 10^^-^^=''° °^ cm^^ for LLSl 134a. The 0.02 dex error 
on the continuum placement is in this case negligible compared to the statistical uncertainty. In 
regards to the metallicity of this gas, we further emphasize that log A^hi> 17.95 is preferred by 
the analysis of the Lyman continuum opacity, despite the symmetric statistical error. If anything, 
metallicity limits derived from this A^hi value may be considered conservative. 

In addition to the H I Lyman series, LLSl 134a exhibits substantial absorption for the D I 
Lyo! and Ly/3 transitions at the expected —82 km offset from the H I lines. Due to blend- 
ing with absorption lines from lower redshift gas, deuterium is not detected in the remaining 
Lyman series lines. Although both of the D I transitions are partially blended with (presumed) 
H I absorption, a simultaneous fit to the hydrogen and deuterium absorption lines (Fig. S2) 
constrains log A^di = 13.26 ± 0.04 and 6di = 10.2 ± 0.8. The latter parameter is fully consis- 
tent with the 6hi value derived from the hydrogen absorption. The deuterium estimate is more 
sensitive to the continuum placement error than hydrogen, particularly at the position of Ly-/3. 
Therefore, we include an additional uncertainty of ~ 0.08 dex, computed by repeating multiple 
times this analysis after having varied the continuum level by ±1 — 3%. When we express 
the Doppler parameter as a function of thermal broadening fethcr and turbulent component feturb, 
the expected ratio of the hydrogen and deuterium b parameters is 6di/&hi = VI ~ 0-50 with 
^ = ^Hi,ther/(^Hi,ther + ^turb) ^^'^ < (f) < 1. The mcasurcd 6di/^hi = 0.66 is within 10% 
of the expected value for turbulent broadening (0 ~ 1). The best-fit D/H ratio for LLSl 134a 
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is log(D/H) = —4.69 ± 0.06, in agreement with previous determinations (Table S3). This 
consistency provides additional confidence on the model for the hydrogen absorption. Due to 
the partial blending of the D I lines with intervening hydrogen, the deuterium column density 
is affected by the additional uncertainty on the line parameters for the hydrogen absorption. 
Including this error, our best estimate becomes log(D/H) = —4.69 ± 0.09. Finally, including 
the uncertainty on the continuum determination, we find log(D/H) = —4.69 ± 0.13. 

2.2 The Lyman Limit System atz = 3.096221 toward Q0956 + 122 

A partial LLS (LLS0956A) was previously identified in low-resolution spectra of Q0956 + 122 
and associated with corresponding metal-line absorption at z — 3.2228 (38). We confirm this 
identification (Fig. 3). An LLS at lower redshift was also reported by Steidel et al. (38), who 
proposed it was associated with a metal-line system at z = 3.1142. Our analysis of the HIRES 
spectra, however, only reveals substantial high-order Lyman series absorption for a system at 
z ~ 3.096 and we adopt this as the redshift for LLS0956B. 

Non-zero flux is detected blueward of the Lyman limit for each system, implying A^hi < 
lO^'^-^ cm-2 for both LLS0956A and LLS0956B. At this optical depth (tll ~ 1), the total A^hi 
value is well established by the Lyman limit opacity and the analysis of unsaturated Lyman 
series lines. 

The hydrogen Lyman series of LLS0956A is detected through to Lyman-21 in our HIRES 
data, together with C IV and Si IV metal lines. By modeling the unsaturated hydrogen tran- 
sitions from Lyman-16 to Lyman-21, we derive the redshift z = 3.223194 ± 0.000002, the 
column density log A^hi= 17.37 ± 0.01, and the Doppler parameter 6hi = 20.4 ± 0.2. Simi- 
larly, for LLS0956B, the hydrogen Lyman series is clearly visible through to Lyman- 15 and, at 
lower signal-to-noise, through to Lyman-21. No other metal lines are detected at the hydrogen 
position, while deuterium is blended in the strongest transitions. Combining both high-order 
saturated and non-saturated H I absorption lines, we model this LLS with a single component at 
z = 3.096221 ±0.000009, with 6hi = 20.2 ±0.8 km s^^ andlogiVHi= 17.18 ±0.03 (Fig. 3 and 
S3). After accounting for an additional 0.02 dex of uncertainty on the continuum placement, the 
best fit column density for LLS0956B becomes log A^hi= 17.18 ± 0.04. The Doppler parameter 
implies a temperature T < 2.5 x 10^ K, typical for photoionized gas. 

3 Metal Line Analysis 

We now discuss upper limits to the column densities of atoms and ions for heavy elements in 
these metal free LLSs. 
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3.1 LLS1134 



LLS1134a is conspicuously free of metal-line absorption, as shown in Fig. 2 where we mark 
for selected ionic transitions the position corresponding to this LLS and to the weaker hy- 
drogen component (z — 3.410883 and z ~ 3.41167 respectively). For both LLS1134a and 
LLS 1134b, we can place upper limits to the ionic column densities using the apparent optical 
depth method (39), with the key limiting factors being the S/N of the spectrum and contamina- 
tion from the Lyct forest (i.e. Lya lines at unrelated wavelengths). For this analysis, we choose 
a velocity window of ±15 km that is wide enough to encompass the expected width for 
metal lines given the hydrogen and deuterium absorption properties. A larger velocity window 
of ±20 km would yield slightly higher (~ 0.06 dex) limits. A summary of the metal column 
densities is presented in Table S2. 

For most of the listed ionic transitions, we place robust upper limits based solely on the vari- 
ance and the uncertainty associated with a 1 — 3% (S/N dependent) error in continuum place- 
ment. Transitions blueward of 5450A, however, lie within the Lya forest which complicates 
the estimates. Of particular interest to the metallicity of LLS 1 134a are the strong transitions of 
Si III 1206 and C III 977. Fig. S4 reveals that there is substantial absorption at the predicted 
wavelength for each transition. This absorption is broad, however, and has peak optical depth 
that is substantially offset from the predicted line-center. It is certain, therefore, that these tran- 
sitions are blended with coincident absorption by the IGM. In this case, an upper limit based on 
the apparent optical depth would overestimate the column densities. 

Better limits to the Si III and C III column densities can be derived using models for these 
absorption lines. We fix the Doppler parameters for the metal lines using the hydrogen and 
deuterium b values. From the H I and D I lines, we infer a temperature (1.6±0.2) x 10^ K 
and an upper limit on the turbulent velocity 6turb < 2.2 km s~^. Carbon and silicon absorption 
lines that arise in gas with these properties would be characterized by be = 6.9 km and 
6si = 5.3 kms^^. These Doppler parameters are within the range commonly observed for 
metal-line transitions (e.g. (40)). A first limit to the column densities of C III and Si III is 
set by demanding that the peak optical depth of each metal line does not exceed the observed 
absorption. This approach, that basically ignores the fact that the opacity is dominated by 
coincident IGM lines, yields upper limits of log iVsim < 11.85 and log A^ciii < 12.45. Lines of 
column density higher than these would exceed the observed absorption at the corresponding 
transitions, as shown in the top panels of Fig. S4. 

A substantial fraction of the opacity at the Si IH and C m frequencies must however arise 
from the two broad Lya forest lines blueward to these metal transitions. We set therefore limits 
to log Ai'siiii and log A^'cm using a two component model that includes both the metal lines and 
the Lya forest. Following this procedure, we find logAi'sini < 11.40 and logAi'ciii < 12.20. 
Since additional Lya lines at both negative and positive velocities are required to account for 
all the observed absorption, the choice of this minimal scenario (only two components) yields 
a quite conservative limit to the metal column densities. 

Inspecting the strongest transitions of LLS 1 134b, the second main hydrogen component at 
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z — 3.41167, we see that Si III and C in are detected. These transitions are commonly found 

in highly ionized LLSs. Again using the apparent optical depth method, we find for the column 
densities of these ions log A^ciii = 13.03 ± 0.02 and logiVsiiii = 12.15 ± 0.02. While a strong 
line is clearly visible at the position of S III for this component, the implied column density 
would greatly exceed that for Si m and C m indicating that the absorption is dominated by 
contamination from the Lya forest. 

3.2 LLS0956B 

Similarly to LLSl 134a, LLS0956B is free of metal line absorption (Fig. 3). Using the apparent 
optical depth method, we determine upper limits to the ionic column densities in the velocity 
interval ±15 km s~^, consistent with the expected metal line width. The metal column den- 
sities for the strongest transitions, including the variance and uncertainty related to continuum 
placement, are listed in Table S2. 

C m A977 is blended in the Lya forest. Supported by the lack of any silicon absorption 
(Si n. Si m, and Si IV) and the absence of C n and C IV, we conclude that carbon is not present 
in this system to limits comparable to the ones inferred from silicon. Therefore, we do not 
include C III A977 in our estimates of the metallicity, but we rely on the Si III upper limit which 
we measure directly with the apparent optical depth method. We note that a larger velocity 
window of ±20 km would result in a slightly higher (0.06 dex) column density limit. We 
also note that the strongest Si n line at A1260 is overwhelmingly contaminated by intervening 
H I absorption and we use the Si n A1526 line to set a limit for this ion. 

4 Metallicity Limits 

Having established the H I column densities and limits to the ionic column densities of Si, C, 
and O, we may set upper limits to the heavy element abundances in the two LLSs. This requires, 
however, a careful consideration of the ionization state of the gas. Generally, the ionization state 
is assessed through a comparison of observed ionic column densities with ionization models. 
Absent the detection of any metal transitions, we must consider physically-motivated scenarios 
for the properties of the gas and thereby estimate limits to the metallicities. 

In a perfectly neutral medium, all of the elements would be in their atomic state (ignoring 
molecular formation) and one could estimate the metallicity directly from the observed column 
densities of the atomic transitions (e.g. H I, C I, O I). In this extreme model one would recover 
log(0/H) -log(O/H)0 = [0/H] < -1.9(-1.4) and [C/H] < -2.1(-1.3) forLLS1134a 
(LLS0956B), which are the most conservative limits possible. In essentially all astrophysical 
environments, however, the gas is irradiated by local and external sources which photoionizes 
at least the outer layers of the cloud. Photons with energy hv <\ Ryd will photoionize all ele- 
ments whose first ionization potential is below that of hydrogen (e.g. C, Si) and harder photons 
may ionized a majority of the H I gas and place the heavy elements in higher ionization states. 
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CoUisional processes, predominantly with electrons and hydrogen nuclei, may further ionize 
the medium for gas at T > 10"^ K. 

Focusing on photoionization first, the population of galaxies and quasars at 2; ~ 3 is known 
to generate an extragalactic UV background (EUVB) which ionizes the majority of baryons in 
the universe (e.g. (41,42)). For gas with the relatively low H I column densities of our LLSs, one 
predicts that the EUVB has substantially ionized the hydrogen gas and any heavy elements that 
are present. To quantitatively assess this process and its impact on the metallicity estimates, we 
have calculated a series of photoionization models using the CLOUDY software package (43). 
We assumed that each LLS may be modeled as a series of layers with constant hydrogen density 
nn, constant and low metallicity, and a plane-parallel geometry. We irradiated this gas with the 
Haardt and Madau EUVB modelj (quasars+galaxies) (14) evaluated at the redshift of each LLS 
(e.g. Ji, = 2.45 X 10^^^ ergs"^ Hz^^ cm~^ at 1 Ryd for z = 3.4) and considered gas with a 
range of rin values. It is standard practice to describe a given photoionization model by the 
ionization parameter U = ^/(nuc) with $ the flux of ionizing photons. Each model was 
restricted to have a total A^hi column density consistent with the observations. 

In Fig. S5, we present the upper limits to the metallicity of the two LLSs as a function 
of nil values. The figure indicates the limit imposed for the range of ions constrained by our 
observations. As a function of nu, each curve shows the metallicity limit from a given ion 
according to the measured upper limit on its column density. The black symbols then trace the 
lowest metallicity imposed by the full set of measurements at each nn value. In all cases we 
compare to the solar abundances of (44). 

At high densities (low U values), the tightest limits are given by the lowest ionization stages 
of the heavy elements (Si"*", 0°, C"*"). These set upper limits to the metallicity of < IO^^Zq. 
We emphasize, however, that such high densities are improbable. A gas with nu > 10^^ cm~^ 
is predicted to have a neutral fraction xhi = nm/nn > 0.1, but all of the LLSs that have 
been studied to date exhibit substantial absorption from higher ionization states of Si, C, and O. 
When analyzed as photoionized gas, one derives ionization parameters log U = —3 to —1 (Fig. 
S6). At the respective redshift for the two LLSs, this implies a volume density nu < 10^^ cm~'^ 
which corresponds to an overdensity relative to the mean baryon density of p/pmcan < 850. 
As nu decreases (and U increases), the heavy elements shift to higher ionization states and 
the limits on the gas metallicity become much more stringent. This is partly a result of tighter 
limits to the ionic column densities but it is primarily because our models are forced to match 
the observed A^hi value. At higher ionization parameter, the gas is more highly ionized and the 
total implied hydrogen column density is correspondingly larger (A^h = Nm/xiu)- Therefore, 
at fixed metallicity this implies higher total column densities for the heavy elements. For the 
ionization parameters typically measured for LLSs (U > 10^'^), we set an upper limit to the 
gas metallicity (based on Si) of Z < 10-^-^{10-^-^)Zq for LLS1134a (LLS0956B). This is a 
conservative value and also the most physically-motivated limit to the metallicity. 

If the intensity of the ionization field is much higher than the assumed EUVB, the lack 

^The results would be slightly different with a different shape for the radiation field, but the differences are 
much smaller than those related to varying the ionization parameter 
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of substantial Si rv and C rv absorption imposes a tight limit to the metallicity, as evident 
from Fig. S5. Extreme scenarios in which most of the metals are in even higher ionization 
states are highly implausible since at log 3> 1 a neutral fraction xm <^ —6 would imply a 
total hydrogen column density A^h 3> 10^"^ cm~^. Further, if these two LLSs lie in proximity 
to a star forming galaxy, the ionization field is dominated by a softer spectrum. To explore 
this possibility, we have calculated a second set of cloudy models adding to the EUVB the 
contribution from a galaxy with star formation rate 100 Mq yr~^. We have further assumed 
that the escape fraction of ionizing radiation is /esc = 0.1 and that the LLSs are at 100 kpc 
from the star forming disks. Under these conditions, the ionization states of the metals do not 
differ substantially from our previous calculation. Indeed, at A^hi = 10^^ cm~^, the inferred 
metallicity for LLS1134a is Z < 10~^ °. Therefore, we conclude that the assumed limits of 
~ 10~^ for LLS1134a and LLS0956B are representative of the underlying gas metallicity in 
different plausible scenarios. 

For gas with T > 2 x 10"^ K, coUisional ionization is an important process for establishing 
the ionization state of the gas. If one assumes coUisional ionization equilibrium, it is straightfor- 
ward to set limits on the gas as a function of the gas temperature (45). For T > 2.5 x 10^ K, the 
limits are yet lower than those suggested by the photoionization modeling. At T ~ 5 x 10"^ K, 
the lack of substantial Si III and C III constrains the metallicity at Z < 10^^{10^^')Zq for 
LLSl 134a (LLS0956B), while at higher temperatures, T ~ 10^ K, the upper limits on the Si IV 
and C IV column densities correspond to Z < 3 x 10-^(10 "5)^0 for LLSl 134a (LLS0956B). 
At temperatures below T ~ 2.5 x 10^ K, we expect both coUisional ionization and photoion- 
ization to contribute, with the latter dominating. 

5 Comparison with Other Observations 

To compare the limits on the metallicity for LLSl 134a and LLS0956B with other observations, 
we search the literature for metallicity estimates of cosmic structures in the redshift interval 
2 < z <5. The metal abundances in solar units as a function of the gas density are summarized 
in Fig. 1. Since we cannot guarantee completeness in our search, this figure offers only a 
qualitative assessment on the interval of metallicity reported in one or multiple studies which 
include large samples. Here we provide a description of the different data shown in the figure. 

The range of densities and metals for the Lya forest is taken from (8), who compile O VI 
absorption statistics at z ~ 2.5 using the line-fitting technique. These metallicities are consistent 
with other studies that are based on different analysis of the IGM (e.g. the pixel optical depth 
method (46)). Upper limits and statistical fluctuations around or below Z ~ 10~^ have been 
also reported at densities p/pmcan ~ 1- These are represented in Fig. 1 using a downward arrow. 

For the LLSs (here loosely defined in the interval log A^hi ~ 10^^ — 10^°'^ cm~^), we search 
in the literature for published metallicity and densities between z — 1.5 — A.S (12, 13, 16, 25, 
47-57). The range plotted in Fig. 1 brackets these observations. For this class of objects we 
also present the logarithmic mean and standard deviation based on our (nearly) complete com- 
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pilation. In addition to the metallicity, we gathered infonnation on the ionization parameters, 
used to produce Fig. S6. ISM metallicity in damped Lya systems has been measured with both 
quasars and 7— ray bursts (GRBs). In both cases, we adopt a number density 0.05 — 10 cm^'^ 
which is consistent with the densities for the ISM atomic phase (58, 59). For the quasar sample, 
we report the range of metallicity between z — 2 — 5 from (5), extended at the lower end with 
the metal poor sample of (60). For the GRB sample we adopt instead the metallicities presented 
in (61). 

Spectroscopy of high redshift galaxies provide abundances from unresolved H II regions 
in star forming galaxies. Here, we adopt the metallicity in two samples of 2; ~ 2 and z ^ 3 
galaxies from (6, 62). We adopt a range of densities typical of nearby H II regions (63). Several 
studies have revealed super-solar metallicity in the broad line regions of quasars. Here, we 
adopt a compilation of metallicity from a large sample of z > 3.7 quasars and we choose a 
density interval typically measured in the broad line region (63). 

Finally, we present the metallicity in Milky Way halo stars. Stellar metallicity is often 
quoted using iron abundance, but the density of heavy elements is dominated by O, C, Ne, and 
other cK-elements. We therefore consider the oxygen abundances presented in (64), extended at 
low metallicity with the star reported by (9). Although the inferred metallicity at the lower end 
suffers from considerable observational uncertainties (64), values at Z ~ 10~^Zq are common 
in large compilations (65) even for stars with extremely low iron abundances. In fact, only a 
single star has been reported at Z ~ 10"^^© (9). 
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Figure SI: LRIS spectrum of Jl 134+5742 at the Lyman limit. Superimposed to the data, we show 
the flux decrement for a A^hi = 10^^"^ cm~^ LLS (blue line) together with the error on the flux (red 
dotted line). In the inset, we show the data at full resolution (gray histogram) and in bins of 50 A 
(crosses). A model of the Lyman limit opacity for log A^hi> 17.90 (green hne) is also shown. 
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Figure S2: Gallery of the Lyman series transitions (black histograms) and best fit models (red lines) 
for LLS1134. The yellow shaded regions represent the 2a errors on the line parameters. LLS1134a and 
LLS1134b are marked by vertical lines, while the dotted black lines indicate the la error on the flux. 
Individual components included in the model are shown with gray lines. 
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Figure S3: Same as Figure S2, but for LLS0956B. 
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Figure S4: Top panels: C III 977 and Si III 1206 transitions for LLS 1134a. Together with the 
HIRES data (black histograms), we display a model for the strongest lines that could be hidden by the 
IGM absorption (red dashed lines). In the bottom panels, we present a single component model (gray 
dashed line) for the IGM absorption that contributes to the opacity at the C III and Si III frequencies. 
Models for these transitions with the adopted column density limits are also shown (gray dotted line). In 
red, we display the combined two-component models, together with the corresponding 2a errors on the 
line parameters. In all cases, we assume a Doppler parameter of 6.9 km s^^ for carbon and 5.3 km 
for silicon, as inferred from the hydrogen and deuterium lines. 
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Figure S5: The two panels describe the metaUicity limits to (top) LLS 1 134a and (bottom) LLS0956B 
as a function of the gas density nn, assuming the Haardt and Madau EUVB radiation field (14). Each 
curve shows the metaUicity limit for a given ion according to the upper limit on its column density. Ele- 
ments C, O, and Si are traced by dotted, dashed, and dash-dot lines respectively. The black symbols then 
trace the lowest metaUicity imposed by the full set of measurements at each nn- The solid gray curves 
in the figure indicate the logarithm of the neutral fraction of the gas xm = nui/nu- 
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Figure S6: Distribution of the ionization parameter U for all the z > 1.5 LLSs reported to date in 
the Uterature. The value U = 10~^ used in this analysis is marked with a dashed Une. 
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Ion 




5v 


Doppler parameter 


Column density 






(km s~^) 


(km s~^) 


(log AT) 


DI 


3.410883 ± 0.000004 





10.2 ±0.8 


13.26 ±0.04 


HI 


3.40997 ±0.00002 


-62 


9.9 ± 1.9 


12.99 ±0.07 


HI 


2.7209 ± 0.0001 




36.9 ± 9.2 


13.05 ±0.12 


HI 


3.410883 ± 0.000004 





15.4 ±0.3 


17.94 ±0.05 


HI 


3.41167 ±0.00001 


54 


18.0 ±0.9 


16.71 ± 0.02 



Table SI: Summary of the best-fit model for the LLS1134 hydrogen and deuterium absorption. For 
each ion we list: the redshift of the absorption, the velocity offset relative to the main hydrogen compo- 
nent, the Doppler parameter and the column density. 
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LLS1134a LLS0956B 



Ion 


A 


log/ 


logiV 




logiV 


1^0 " 




(A) 




(cm ) 


(mA) 


(cm ) 


(mA) 


CII 


1334.5323 


-0.8935 


< 12.45 


< 5.7 


< 12.26 


< 3.7 


cm 


977.0200 


-0.1180 


< 12.20 


< 10.2 






CIV 


1548.1950 


-0.7194 


< 12.09 


< 5.0 


< 12.30 


< 8.1 


01 


1302.1685 


-1.3110 


< 12.76 


< 4.2 


< 12.50 


< 2.3 


Si II 


1260.4221 


0.0030 


< 11.30 


< 2.8 






Sin 


1526.7066 


-0.8962 


< 12.27 


< 4.9 


< 12.45 


< 7.4 


Sim 


1206.5000 


0.2201 


< 11.40 


< 5.4 


< 11.21 


< 3.5 


Si IV 


1393.7550 


-0.2774 


< 11.69 


< 4.4 


< 11.83 


< 6.1 


Fen 


1608.4511 


-1.2366 


< 12.59 


< 5.2 


< 12.72 


< 7.0 



" The rest-frame equivalent width (Wq) is computed for the hnear portion of the curve of growth. 



Table S2: 2a upper limits on the column density of metal ions and rest-frame equivalent widths. Rest 
frame wavelengths and oscillator strengths (/) from (66) are also listed. 
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Quusur 




log u / ri 


log 








_A fin + n 04 


1 q AO -1- n m 




Q0913 + 072 


2.61843 


-4.56 ±0.04 


20.34 ±0.04 


-2.37" 


Q1009 + 299 


2.50357 


-4.40 ± 0.07 


17.39 ± 0.06 


-2.5* 


Jl 134+5742 


3.41088 


-4.69 ±0.13 


17.95 ± 0.05 


< -4.2*'^ 


Q1243 + 3047 


2.52566 


-4.62 ± 0.05 


19.73 ±0.04 


-2.79" 


SDSSJ1558 - 0031 


2.70262 


-4.48 ± 0.06 


20.67 ±0.05 


-1.47" 


Q1937 1009 


3.57220 


-4.48 ± 0.04 


17.86 ±0.02 


-2.7,-1.9*'"' 


Q2206 - 199 


2.07624 


-4.78 ±0.09 


20.43 ± 0.04 


-2.04" 



" Metallicity obtained using oxygen. ^ Metallicity obtained using silicon. This work. The authors report metal- 



licities for two components. Data from (2, 5, 16, 49-51, 67, 68). 



Table S3: Deuterium abundances in z > 2 absorption line systems. 
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